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ABSTRACT 


\ 

V 

Electric  field  spectrum  measurements  from  the  Plasma  Wave  Instru¬ 
ment  on  the  Dynamics  Explorer-1  spacecraft  are  used  to  study  the  local 
electron  density  at  high  altitudes  in  the  northern  polar  cap  region. 

The  electron  density  is  determined  from  the  upper  cutoff  of  whistler 
mode  radiation  at  the  electron  plasma  frequency.  This  technique  typi¬ 
cally  gives  the  electron  density  with  an  accuracy  better  than  12%.  The 
median  electron  density  over  the  polar  cap  at  L  greater  than  10  is'^X^ 
found  to  vary  from  35.2  ±  8.5  cm-3  at  2.1  Rg  to  0.99  ±  0.51  cm“3  at 
4.66  Rg.  Variations  up  to  a  factor  of  four  occur  from  these  median 
values . 

The  steady  state  radial  outflow  model  is  examined  for  consistency 
with  the  observed  density  profile.  The  model  predicts  an  inverse 
dependence  on  the  product  of  the  plasma  velocity  and  the  cube  of  the 
radial  distance.  A  power  law  fit  to  the  radial  variation  of  the  elec¬ 
tron  density  yields  an  exponent  of  -3.85  ±  0.32,  which  for  the  radial 
outflow  model  implies  a  flow  velocity  increasing  nearly  linearly  with 
Increasing  radial  distance.  Comparison  of  the  observed  electron  den¬ 
sities  with  theoretical  polar  wind  densities  yields  consistent  results 
up  to  2.8  Rg,  although  the  steady-state  solutions  for  supersonic  H+ 
outflow  are  slightly  lower  than  the  observed  median  densities.  Recent 
DE-I  measurements  contradict  the  polar  wind  predictions  of  ion 


composition  and  energies,  although  a  thermal  polar  wind  component  in 
the  polar  cap  plasma  is  also  present. 


A  comparison  of  observed  electron  densities  with  low-altitude  den¬ 
sity  profiles  from  the  Alouette  11  and  Isis-1  spacecrafts  Indicates 
that  a  transition  in  the  slope  of  the  electron  density  profile  occurs 
at  1.16  Rg  and  that  another  transition  occurs  between  1.55  Rg  and  2.0 
Rg.  These  transitions  are  probably  due  to  changes  in  the  basic  pro¬ 
cesses  of  radial  plasma  transport  in  the  polar  cap  region.  At  the 
first  transition  the  plasma  transport  is  believed  to  change  from  a 
combined  subsonic  flow  and  diffusive  equilibrium  distribution  below 
1.16  Rg  to  a  collision-dominated  supersonic  flow  above  1.16  Rg.  At  the 
second  transition  the  plasma  transport  is  believed  to  change  from  a 
collision-dominated  outflow  below  1.5  Rg  to  a  collisionless  flow  above 
2.0  Rg. 
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I.  INTRODUCTION 


In  this  report  the  electron  density  profile  as  a  function  of 
radial  distance  is  determined  using  electric  field  spectrum  measure¬ 
ments  from  the  Plasma  Wave  Instrument  on  the  Dynamics  Explorer- 1  space¬ 
craft.  The  region  of  interest  in  this  study  is  the  earth's  northern 
polar  cap  region  where  geomagnetic  field  lines  open  into  the  mag¬ 
netotail  and  a  steady  outward  plasma  flow  occurs.  The  polar  cap  region 
is  defined  to  be  the  region  of  the  magnetosphere  poleward  of  the  auro¬ 
ral  field  lines.  For  purposes  of  this  study,  the  northern  polar  cap 
region,  was  taken  to  be  that  region  of  the  magnetosphere  with  L-shell 
values  greater  than  10,  corresponding  to  invariant  latitudes  greater 
than  71.6°  (A  ■  cos”1  i/l/L) . 

The  Dynamics  Explorer-1  spacecraft,  launched  on  August  3,  1981,  is 
in  an  eccentric  polar  orbit  with  perigee  and  apogee  geocentric  radial 
distances  of  1.09  %  and  4.66  Rg,  respectively.  Because  of  the  eccen¬ 
tricity  and  the  latitudinal  precession  of  the  orbit,  this  spacecraft 
provides  electric  and  magnetic  field  measurements  over  a  region  of  the 
polar  magnetosphere  not  previously  studied.  Figure  1  illustrates  the 
orientation  of  the  DE-1  orbit  in  September,  1981  and  the  precession  of 
the  apogee  over  the  five-month  period  involved  in  this  study,  from 
September,  1981  to  February,  1982.  The  latitudinal  precession  of  the 
DE-1  apogee  is  108s  per  year.  This  precession  is  responsible  both  for 


the  variations  in  the  altitude  of  the  spacecraft  in  the  polar  cap 
region  and  for  the  variations  in  the  length  of  time  the  spacecraft 
spends  in  this  region.  During  the  five-month  period  covered  in  this 
study,  the  geocentric  radial  distance  of  DE-1  in  the  northern  polar  cap 
region  varied  from  2.0  Rg  to  4.66  Rg,  or  6378  km  to  23,300  ka  above  the 
earth's  surface.  This  altitude  range  is  far  below  the  high  altitudes 
surveyed  by  the  Hawkeye  spacecraft  and  above  the  low  altitudes  surveyed 
by  the  Alouette-lsis  spacecraft  and  the  S3-3  spacecraft. 

Prior  to  the  1960's,  electron  density  measurements  up  to  the 
height  of  the  F2  peak  were  obtained  with  the  use  of  ground-based 
ionosondes  which  recorded  the  reflection  of  radio  wave  pulses  from  the 
ionosphere  [Calvert,  1966].  To  study  the  electron  density  profile 
above  the  F2  peak,  topside  sounders,  similar  in  concept  and  operation 
to  ground-based  ionosondes,  were  designed  for  operation  on  orbiting 
satellites.  Four  topside  sounder  satellites  of  the  Alouette-lsis 
series  were  launched  in  the  decade  following  1961  and  provided  electron 
density  measurements  from  the  F2  peak  up  to  about  3300  km  [Jackson  and 
Warren,  1969;  Florida,  1969;  Coordinated  Ionospheric  and  Magnetospheric 
Observations  from  the  Isis  2  Satellite,  1980].  These  topside  sounder 
satellites  were  able  to  obtain  electron  density  profiles  over  the 
entire  polar  cap  region  at  altitudes  below  3500  km. 

Later,  two  spacecraft,  S3-3  and  Hawkeye  1,  were  launched  into  even 
higher  altitude  orbits  over  the  polar  cap.  The  S3-3  satellite, 
launched  in  the  summer  of  1976,  carried  a  Langmuir  probe  Instrument  to 
measure  electron  number  densities  and  density  fluctuations  along  the 


satellite  orbit  [Mozer  et  al.,  1979].  The  S3-3  spacecraft  was  placed 
in  an  elliptical  polar  orbit  with  an  apogee  of  8000  km.  However, 
although  the  data  are  available,  no  electron  density  measurements  in 
the  polar  cap  region  have  yet  been  published.  The  Hawkeye  I  space¬ 
craft,  launched  on  June  3,  1974,  was  placed  in  a  highly  eccentric  orbit 
with  a  polar  apogee  of  20.5  Rg.  Due  to  the  precession  of  its  apogee 
during  its  four-year  lifetime,  the  Hawkeye  spacecraft  was  able  to 
obtain  electric  and  magnetic  field  measurements  and  corresponding 
charge  particle  measurements  at  large  radial  distances  up  to  20  Rg  in 
the  polar  cap  region  [Gurnett  and  Frank,  1978].  However,  because  of 
the  eccentricity  of  its  orbit,  there  was  a  gap  in  the  Hawkeye  data 
coverage  over  the  polar  cap  below  approximately  8  Rg.  Although  density 
measurements  at  higher  altitudes  in  the  polar  cap  region  could  be 
derived  from  the  Hawkeye  electric  field  measurements  [Calvert,  1981], 
to  date  no  such  results  have  been  published. 

The  Plasma  Have  Instrument  (PHI)  on  Dynamics  Explorer-1  has  been 
described  in  detail  in  a  previous  report  [see  Shawhan  et  al.,  1981]. 

The  instrument  is  designed  to  measure  AC  electric  fields  over  a  fre¬ 
quency  range  of  1  Hz  to  2  MHz  and  an  amplitude  range  of  0.03  pV  m~ 1  to 
100  mV  m-1.  AC  magnetic  fields  and  quasi-static  electric  fields  are 
also  measured.  The  antenna  system  used  to  detect  electric  fields  con¬ 
sists  of  a  pair  of  electric  antennas  mounted  orthogonally  to  each 
other.  One  of  these  antennas  is  9  m  long  and  is  oriented  parallel  to 
the  spin  axis  and  the  other  is  200  m  long  and  is  oriented  perpendicular 
to  the  spin  axis.  Signals  from  these  sensors  are  processed  by  two  of 
the  PHI  receiver  systems,  the  Step  Frequency  Correlator  (SFC)  and  the 


Low  Frequency  Correlator  (LFC).  The  Step  Frequency  Correlator  consists 
of  a  correlator  and  two  identical  Step  Frequency  Receivers  (SFR)  with 
four  frequency  bands,  each  containing  32  frequency  steps.  The  SFC 
scans  a  frequency  range  of  100  Hz  to  400  kHz  once  every  32  seconds. 

The  Low  Frequency  Correlator  consists  of  a  correlator  and  two  identical 
spectrum  analyzers,  each  with  eight  filters.  The  LFC  scans  a  frequency 
range  of  1.78  Hz  to  100  Hz  once  every  32  seconds.  The  resulting  elec¬ 
tric  field  amplitude  measurements  are  displayed  in  frequency-time  spec¬ 
trogram  form.  Each  spectrogram  covers  a  two-hour  time  interval  and 
spans  a  frequency  range  of  1.78  Hz  to  400  kHz. 


II.  METHOD  OF  ANALYSIS 


A.  Plasma  Wave  Modes  in  the  Polar  Cap 

Electric  field  spectrograms  obtained  from  the  electric  field 

amplitude  measurements  of  the  PWI  illustrate  features  of  the  various 

modes  of  plasma  waves  observed  by  DE-1  for  frequencies  in  the  vicinity 

of  the  electron  gyrofrequency  and  the  electron  plasma  frequency.  Four 

plasma  wave  modes  are. predicted  by  cold  plasma  theory  in  this  frequency 

range  [ Stix,  1962].  These  modes  are  the  free  space  L-0  (left-hand 

polarized,  ordinary)  mode,  the  free  space  R-X  (right-hand  polarized, 

extraordinary)  mode,  the  Z-mode  (extraordinary)  and  the  whistler  mode 

(right-hand  polarized).  The  free  space  L-0  mode  has  a  low  frequency 

cutoff  at  the  electron  plasma  frequency.  The  free  space  R-X  mode  has  a 

fB  fff  2  2  W2 

low  frequency  cutoff  at  f(R-O)-  -g®-  +  ((-g®-)  +  fp)  .  The  Z-mode  is 

bounded  at  low  frequencies  by  f(L-O)  -  -  ■=&  +  ((-A 2  +  f2)1/2  and  at 

2  2  ^2  21/2 
high  frequencies  by  the  upper  hybrid  resonance,  fmm  ■  (fg  +  fp)  . 

The  whistler  mode  is  confined  to  frequencies  below  the  electron  gyro- 
frequency  or  the  electron  plasma  frequency,  whichever  is  smaller. 

These  four  plasma  wave  modes  have  been  associated  with  electro¬ 
magnetic  plasma  wave  emissions  found  in  the  auroral  zone  and  polar  cap 
region  [Gurnett  et  al.,  1983].  Broadband  electromagnetic  emissions 
propagating  below  the  electron  gyrofrequency  in  the  auroral  zone  (where 
fg  <  fp)  are  called  aur  --1  his  These  emissions  are  believed  to  be 
propagating  in  the  whistler  mode  [Gurnett  et  al.,  1983]  because  this  is 


the  only  electromagnetic  plasma  wave  mode  propagating  in  this  frequency 
range.  Auroral  hiss  is  generated  by  low  energy  (100  eV  to  10  keV) 
electrons  moving  upward  along  auroral  field  lines  from  a  region  of 
electron  acceleration  between  1.7  Rg  and  1.9  Rg  [Lin  et  al.t  1983]. 
Intense,  high-frequency  emissions  called  auroral  kilometric  radiation 
(AKR)  have  also  been  identified  on  the  PWI  electric  field  spectrograms. 
This  radiation  is  generated  along  auroral  field  lines  between  2  Rg  and 
4  Rg  [Gurnett,  1974;  Gallagher  and  Gumett,  1979]  and  escapes  from  the 
Earth's  magnetosphere.  Recent  polarization  measurements  of  AKR  by 
Shawhan  and  Gurnett  [ 1982]  have  determined  that  the  dominant  component 
of  AKR  is  right-hand  polarized  and  the  radiation  propagates  in  the  free 
space  R-X  mode.  Broadband  Z-mode  radiation  has  been  identified  in  the 
auroral  zone  by  Calvert  [1981]  and  Gurnett  et  al.  [1983]. 


B.  A  Typical  Spectrogram:  Day  309,  1981 
A  typical  example  of  an  electric  field  spectrogram  illustrating 
these  various  wave  modes  is  shown  in  Figure  2.  The  trajectory  of  the 
spacecraft  in  a  magnetic  meridian  plane  is  shown  in  the  lower  left 
corner  of  the  spectrogram.  In  this  case,  the  spacecraft  is  approaching 
the  nightslde  auroral  region  from  a  radial  distance  of  4.2  Rg  and  79° 
magnetic  latitude  and  entering  the  auroral  zone  at  34°  magnetic  lati¬ 
tude  and  a  radial  distance  of  2.8  Rg.  The  intense  emissions  in  the 
high  frequency  range  are  auroral  kilometric  radiation.  The  apparent 
lower  frequency  cutoff  of  auroral  kilometric  radiation  is  an  instru¬ 
mentation  effect;  the  low  frequency  cutoff  of  the  Step  Frequency 
Receivers'  fourth  frequency  band  is  locate  '  at  58  kHz.  The  emissions 


T 


In  the  IS  -  60  kHz  range  are  Z-mode  radiation,  electromagnetic  emis¬ 
sions  with  a  characteristic  upper  frequency  cutoff  located  at  the  upper 
hybrid  resonance  frequency  [Gurnett  et  al.,  1983].  In  this  spectro¬ 
gram,  the  Z-mode  radiation  exhibits  a  sharply  defined  upper  frequency 
cutoff  at  the  electron  gyrofrequency  on  the  poleward  side  of  the  auro¬ 
ral  zone.  In  this  region,  which  has  a  low  electron  density,  the  elec¬ 
tron  plasma  frequency  is  much  lower  than  the  local  electron  gyrofre¬ 
quency.  Under  these  circumstances,  the  upper  hybrid  resonance  fre¬ 
quency  is  approximately  equal  to  the  local  electron  gyrofrequency  and 
the  latter  becomes  the  upper  frequency  cutoff  for  the  Z-mode. 

The  intense  emissions  below  20  kHz  which  spread  out  from  the  auro¬ 
ral  region  toward  the  pole  are  auroral  hiss  emissions.  Auroral  hiss  is 
commonly  observed  between  10  Hz  and  100  kHz.  Figure  2  illustrates  the 
funnel-shaped  lower  frequency  cutoff  of  auroral  hiss  emissions  in  the 
auroral  zone,  a  feature  that  is  frequently  but  not  always  observed  at 
the  altitude  of  the  DE-1  orbit  [Gurnett  et  al.,  1983].  Auroral  hiss  is 
not  confined  to  the  auroral  zone  but  spreads  out  over  a  wide  range  of 
latitudes  both  poleward  and  equatorward  of  the  auroral  zone.  Auroral 
hiss  emissions  are  also  found  in  the  dayside  polar  cusp  and  sometimes 
exhibit  a  latitudinal  asymmetry  in  the  direction  of  the  pole  (see 
Figure  3),  resulting  in  a  nearly  continuous  band  of  emissions  over  the 
polar  cap  region  (see  Figures  4  and  5). 

Since  auroral  hiss  is  a  whistler  mode  emission,  this  radiation  has 


a  characteristic  upper  frequency  cutoff  at  either  the  electron  gyro¬ 
frequency  or  the  electron  plasma  frequency,  whichever  is  smaller  [Stix, 
1962;  Chap.  2;  Chen,  1974,  Chap.  4],  In  the  polar  cap  region,  the 


electron  plasma  frequency  Is  generally  much  smaller  than  the  electron 
gyrofrequency.  The  electron  plasma  frequency  then  gives  the  upper  fre¬ 
quency  cutoff  for  the  whistler  mode.  In  Figure  2,  the  plasma  frequency 
cutoff  is  well-defined  and  varies  smoothly  from  16  to  19  kHz.  It 
distinguishes  the  auroral  hiss  emissions  below  from  the  Z-mode  radi¬ 
ation  above. 

Spectrograms,  such  as  the  one  in  Figure  2,  are  constructed  from 
the  digitized  data  output  of  the  Step  Frequency  Receivers.  To  trace 
the  plasma  frequency  cutoff,  the  digitized  data  are  displayed  on  a 
video  monitor  in  spectrogram  form.  A  manually  controlled  cursor  is 
then  used  to  trace  the  plasma  frequency  cutoff  on  the  spectrogram  and  a 
frequency  value  is  recorded  for  each  sweep  of  the  Step  Frequency 
Receiver  (every  32  seconds).  In  addition,  orbit  parameters  of  the 
spacecraft,  such  as  magnetic  latitude  and  radial  distance,  are  recorded 
for  each  32-second  time  Interval.  From  a  knowledge  of  the  electron 
plasma  frequency,  the  electron  number  density  can  be  derived,  using  the 
well-known  relation: 

n  -  (f  /9)2 
e  p 

where  ng  is  the  electron  number  density  in  cm-3  and  fp  is  the  electron 
plasma  frequency  in  kHz. 

In  addition,  a  quality  index  is  assigned  to  each  data  point.  The 
quality  index  is  used  to  identify  data  of  questionable  reliability  in 
cases  where  the  Intensity  of  the  whistler  mode  radiation  approaches  the 


intensity  of  the  background  noise  and  the  upper  frequency  cutoff  is  no 
longer  clearly  defined.  The  use  of  the  quality  control  index  makes  it 
possible  to  study  the  effect  of  marginally  reliable  data  on  the  density 
profile  and  to  locate  spurious  data  points. 

For  spectrograms  in  which  the  upper  frequency  cutoff  is  not 
clearly  defined,  contrast  enhancement  techniques  are  employed  to  accent 
the  plasma  frequency  cutoff  and  minimize  background  noise  and  Z-mode 
radiation.  Contrast  enhancement  techniques  are  helpful  in  locating  the 
frequency  cutoff  when  the  intensity  of  the  Z^mode  radiation  approaches 
the  intensity  of  the  whistler  mode  radiation,  thereby  partially  obscur¬ 
ing  the  plasma  frequency  cutoff.  This  effect  is  evident  after  0550  in 
Figure  2.  A  16-point  sliding  grey  scale  provides  the  capability  of 
enhancing  the  intensity  of  segments  of  the  frequency  range  by  raising 
the  minimum  intensity  levels  above  the  maximum  intensity  of  the  back¬ 
ground  noise  and  lowering  the  maximum  intensity  levels  displayed  on  the 
monitor  screen.  The  remaining  narrower  intensity  range  is  spanned  by 
16  shades  from  white  (lowest  Intensity)  to  black  (highest  intensity). 
Such  contrast  enhancement  techniques  are  necessary  since  an  accurate 
determination  of  the  electron  number  density  depends  on  the  measurement 
of  a  well-defined  plasma  frequency  cutoff. 

The  combined  use  of  contrast  enhancement  techniques  and  tracing 
the  frequency  cutoff  with  a  cursor  permits  the  determination  of  the 
frequency  cutoff  to  a  high  degree  of  accuracy.  The  top  three  SFR  fre¬ 
quency  bands  (900  Hz  -  400  kHz)  are  divided  into  96  frequency  steps. 

The  cursor  can  locate  the  frequency  cutoff  within  one  frequency  step. 

At  10  kHz,  a  typical  frequency  cutoff  value,  the  frequency  cutoff  can 


be  determined  within  ±600  Hz,  corresponding  to  an  uncertainty  in  the 
calculated  number  density  of  ±122.  This  uncertainty  is  negligible 
compared  to  the  observed  spread  in  density  values  due  to  diurnal  and 
latitudinal  variations,  as  well  as  anticipated  large  scale  variations, 
due  to  magnetic  activity  and  solar  effects. 

An  additional  source  of  error  comes  from  the  identification  of  the 
cutoff  frequency  with  the  electron  plasma  frequency.  Hot  plasma 
effects  such  as  cyclotron  damping  and  Landau  damping  can  cause  the 
whistler  mode  wave  to  cut  off  at  frequencies  below  the  electron  plasma 
frequency.  These  effects,  however,  are  expected  to  be  small.  Because 
of  the  small  wave  normal  angle  of  the  whistler  mode  at  resonance  and 
the  subsequently  high  index  of  refraction  required  for  Landau  damping, 
this  effect  will  only  be  significant  in  a  narrow  band  of  frequencies 
Just  below  the  plasma  frequency.  Only  particles  with  resonant  frequen¬ 
cies  in  this  band  will  interact  with  the  wave.  For  frequencies  in  the 
range  of  fp  ±0.2  kHz,  the  Landau  damping  effect  would  imply  an  error 
in  the  calculated  number  density  of  approximately  42  at  fp  -  10  kHz. 
Cyclotron  damping  will  have  a  significant  effect  on  the  cutoff  fre¬ 
quency  only  in  a  region  near  the  electron  gyrofrequency.  A  statistical 
comparison  of  the  electron  gyrofrequency  and  the  observed  cutoff  fre¬ 
quency  on  the  electric  field  spectrograms  indicates  that  the  cutoff 
frequency  was  equal  to  or  greater  than  0.9  fg  for  approximately  22  of 
all  data  points  in  this  study  (see  Section  IV).  Consequently,  the 
cyclotron  damping  effect  implies  an  error  in  the  calculated  number  den¬ 
sity  for  22  of  all  density  data  points.  The  cyclotron  damping  effect 


is  negligible  below  3.5  Re  and  is  expected  to  be  smaller  than  the 
Landau  damping  effect  at  all  radial  distances. 

The  method  of  determining  the  electron  number  density  from  the 
plasma  frequency  cutoff  has  certain  inherent  limitations.  The  plasma 
frequency  can  be  determined  with  the  aforementioned  accuracy  only  for 
those  spectrograms  where  auroral  hiss  is  present  and  where  the  cutoff 
is  reasonably  well-defined.  Of  the  more  than  750  electric  field  spec¬ 
trograms  examined  in  this  study,  only  40%  have  auroral  hiss  with  well- 
defined  upper  frequency  cutoffs  for  a  significant  fraction  of  the  two- 
hour  time  interval  of  the  spectrogram.  The  quality  of  the  upper  fre¬ 
quency  cutoff  on  these  spectrograms  varies  with  the  intensity  of  the 
whistler  mode  radiation  as  well  as  with  the  intensity  of  the  background 
noise  and  the  other  plasma  wave  modes  (particularly  the  Z-mode)  propa¬ 
gating  near  the  upper  frequency  cutoff.  The  quality  of  the  frequency 


III.  REPRESENTATIVE  ELECTRIC  FIELD  SPECTROGRAMS 


To  provide  a  feeling  for  the  type  of  data  used  and  the  diffi¬ 
culties  encountered  five  representative  passes  have  been  selected  for 
discussion.  These  passes  were  selected  to  illustrate  a  range  of  condi¬ 
tions,  varying  from  very  smooth  well-defined  electron  density  profiles 
to  highly  Irregular  profiles  with  large  variations  on  a  short  time 
scale. 


A.  Day  278,  1981 

The  spectrogram  for  this  pass  is  shown  in  Figure  3.  The  electron 
density  profile  is  smooth  and  well-defined  up  until  0840  as  the  satel¬ 
lite  nears  the  pole.  The  spectrogram  is  from  early  October,  1981,  when 
the  satellite's  apogee  had  precessed  about  15s  to  the  dayside  of  the 
pole.  The  panel  in  the  lower  left-hand  corner  of  the  spectrogram  indi¬ 
cates  the  satellite's  trajectory.  DE-1  is  emerging  from  the  dayside 
polar  cusp  at  4.3  Rg  with  a  magnetic  latitude  of  59°  and  approaching 
the  pole  at  4.6  Rg  and  78°  magnetic  latitude  at  0858  where  there  is  a 
data  gap.  Intense  AKR  is  again  evident  above  100  kHz  and  Z-mode  radi¬ 
ation,  less  intense  than  it  was  in  Figure  2,  is  found  between  10  kHz 
and  20  kHz  at  0755.  The  auroral  hiss  emissions  are  latitudlnally 
asymmetric,  spreading  out  over  the  polar  cap  region  in  the  direction  of 
the  pole.  However,  the  funnel-shaped  lower  frequency  cutoff  of  the 
auroral  hiss  in  the  auroral  zone  is  not  evident  in  this  spectrogram. 
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Below  78°  magnetic  latitude,  the  electron  plasma  frequency  Is  smoothly 
varying,  decreasing  with  increasing  radial  distance.  Exceptions  are 
the  sharp  spike  in  the  plasma  frequency  cutoff  at  0752  and  several 
spikes  after  0845,  indicating  a  sudden  increase  in  the  electron  number 
density  by  more  than  a  factor  of  two.  After  about  0840  the  auroral 
hiss  becomes  weak  and  the  upper  cutoff  becomes  irregular  and  difficult 
to  identify. 


B.  Day  047,  1982 

The  two  spectrograms  in  Figures  4  and  5  were  chosen  to  illustrate 
one  continuous  pass  of  DE-1  over  the  polar  cap.  This  pass  occurred  in 
mid-February  of  1982  when  the  spacecraft  apogee  had  precessed  into  the 
nightslde  auroral  zone.  In  Figure  4,  the  spacecraft  is  in  and  emerging 
from  the  night-side  auroral  zone  as  illustrated  by  the  panel  in  the 
bottom  left-hand  corner  of  the  spectrogram.  It  passes  through  apogee 
at  approximately  0350  and  emerges  from  the  auroral  zone  shortly  after 
0405,  with  a  range  in  magnetic  latitude  of  34°  to  69°.  An  abrupt  drop 
in  the  plasma  density  is  evident  in  the  auroral  zone  from  about  0400  to 
0405.  This  density  depletion  is  believed  to  correspond  to  the  auroral 
plasma  cavity  discussed  by  Calvert  [1981].  In  the  polar  cap  region, 
the  plasma  frequency  cutoff  of  the  whistler  mode  radiation  is  sharply 
defined  and  constant  at  approximately  10  kHz.  The  Z-mode  radiation 
above  the  auroral  hiss  is  less  intense  and  again  exhibits  a  sharp 
upper-frequency  cutoff  at  the  electron  gyrof requency. 

The  spectrogram  in  Figure  5  is  a  continuation  of  the  two-hour 
spectrogram  in  Figure  4.  The  spacecraft  is  in  the  northern  polar  cap 
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region  for  most  of  this  time  interval,  passing  over  the  pole  at  approx¬ 
imately  0600.  The  vertical  white  lines  on  the  spectrogram  indicate 
short  data  gaps.  The  radial  distance  of  the  spacecraft  for  this  spec¬ 
trogram  varies  from  4.4  Rg  at  a  magnetic  latitude  of  62°  on  the  night- 
side  of  the  pole  to  1.8  %  at  a  magnetic  latitude  of  44°  on  the  dayside 
of  the  pole.  The  spacecraft  enters  the  polar  cusp  at  approximately 
0640,  as  Indicated  by  the  onset  of  broadband  auroral  hiss  emissions  at 
that  time.  The  plasma  frequency  cutoff  in  Figure  5  is  again  sharply 
defined  and  is  generally  constant  at  about  16  kHz  before  the  spacecraft 
crosses  the  pole.  However,  there  is  an  abrupt  12%  decrease  In  magni¬ 
tude  at  a  magnetic  latitude  of  84°  on  the  dayside  of  the  pole  (at  0604) 
and  a  steady  increase  in  the  plasma  frequency  after  0625  with  decreas¬ 
ing  radial  distance.  Generally,  though,  the  plasma  frequency  in 
Figures  4  and  5  exhibits  few  variations  on  a  short  time  scale  and  the 
resulting  electron  density  profile  for  this  orbital  pass  is  smoothly 
varying. 


C.  Day  046,  1982 

However,  the  plasma  frequency  is  not  smoothly  varying  for  the 
orbital  pass  which  immediately  precedes  the  pass  illustrated  in  Figures 
4  and  5.  In  Figure  6  the  spacecraft  is  in  the  same  region  of  the  polar 
magnetosphere  as  it  is  in  Figures  4  and  5,  only  6  hours  earlier.  As 
indicated  in  the  panel  in  the  lower  right-hand  corner,  DE-1  is  emerging 
from  the  nlghtside  auroral  zone  into  the  northern  polar  cap  region.  In 
this  two-hour  interval,  the  radial  distance  of  the  spacecraft  varies 
from  4.4  Rg  at  a  magnetic  latitude  of  47°  on  the  nlghtside  of  the  pole 
to  1.8  Rg  at  a  magnetic  latitude  of  51°  on  the  dayside  of  the  pole. 
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We  believe  that  the  radiation  seen  below  40  kHz  at  high  magnetic 
latitudes  in  Figure  6  is  electromagnetic.  The  radiation  propagates  at 
frequencies  below  the  electron  plasma  frequency  in  a  region  of  the 
polar  magnetosphere  where  fp  «  fg.  Whistler  mode  radiation  is  the 
only  electromagnetic  plasma  wave  mode  that  can  propagate  in  this  fre¬ 
quency  range.  This  whistler  mode  radiation,  however,  is  not  the  usual 
type  of  auroral  hiss  associated  with  the  auroral  zone  since  it  is 
believed  to  be  generated  at  high  magnetic  latitudes  poleward  of  the 
auroral  field  lines.  This  radiation  may  be  associated  with  the  polar 
cap  aurora  or  other  aurora-like  processes  over  the  polar  cap. 

Since  this  radiation  is  whistler  mode  radiation,  the  upper  fre¬ 
quency  cutoff  visible  in  Figure  6  is  the  electron  plasma  frequency. 
Unlike  the  plasma  frequency  cutoff  in  Figures  4  and  5,  the  plasma  fre¬ 
quency  cutoff  in  Figure  6  is  highly  variable  on  time  scales  of  8 
minutes  to  less  than  one  minute,  indicating  variations  in  the  electron 
number  density  of  an  order  of  magnitude  or  more  on  spatial  scales  of 
50  to  1400  km. 


D.  Day  039,  1982 

The  variability  in  the  character  of  the  electron  density  profile 
between  successive  orbital  passes  of  the  spacecraft  is  again  illus¬ 
trated  in  Figures  7  and  8.  The  spectrograms  in  these  figures  are  from 
successive  passes  of  the  spacecraft,  emerging  from  the  nlghtside  auro¬ 
ral  zone  in  early  February,  1982.  In  Figure  7,  the  spacecraft  exits 
the  auroral  zone  at  0400  at  a  radial  distance  of  4.6  Rg  and  a  magnetic 
latitude  of  approximately  45s  and  passes  through  apogee  in  the  polar 
cap  region  about  10  minutes  later.  In  this  spectrogram,  the  electron 


plasma  frequency  varies  smoothly  between  9  khz  and  13  kHz  with  few 
abrupt  changes  on  a  short-time  scale,  indicating  changes  in  the  elec¬ 
tron  number  density  by  a  factor  of  2  or  less. 

On  the  following  pass  (Figure  8)  the  spacecraft  emerges  from  the 
nightside  auroral  zone  at  1140  (approximately  7.5  hours  later)  at  a 
radial  distance  of  4.5  Rg  and  a  magnetic  latitude  of  50°.  The  radial 
distance  of  the  spacecraft  in  the  polar  cap  region  for  this  pass  varies 
from  4.5  %  at  a  magnetic  latitude  of  50°  to  3  %  at  a  magnetic  lati¬ 
tude  of  78°  on  the  nightside  of  the  pole,  very  similar  to  the  preceding 
pass.  However,  unlike  the  preceding  pass,  the  electron  plasma  fre¬ 
quency  cutoff  is  not  smoothly  varying.  There  is  a  general  tendency  for 
the  plasma  frequency  to  increase  in  magnitude  with  decreasing  radial 
distance  as  the  spacecraft  nears  the  pole,  as  it  did  in  Figure  7.  But 
more  striking  are  the  large  variations  in  the  magnitude  of  the  plasma 
frequency  on  time  scales  of  5  minutes  or  less,  indicating  variations  in 
the  electron  number  density  up  to  one  order  of  magnitude  on  spatial 
scales  of  less  than  700  km.  The  electron  density  profile  for  this  pass 
is  not  as  smoothly  varying  as  it  was  in  Figure  7. 


IV.  THE  ELECTRON  DENSITY  PROFILE 


To  obtain  the  electron  density  profile  in  the  northern  polar  cap 
region,  the  electron  plasma  frequency  was  numerically  determined  from 
the  analysis  of  nearly  one  hundred  2-hour  spectrograms  over  the  polar 
cap,  selected  from  the  five-month  time  interval  from  September,  1981 
to  February,  1982.  This  Interval  was  characterized  by  relatively  high 
magnetic  activity  (kp  =3).  The  electron  number  density  was  calculated 
from  the  plasma  frequency  using  the  relation:  rig  -  (fp/9)2  where  ne  is 
in  cm-3  and  fp  is  in  kHz.  Figure  9  is  a  semilog  plot  of  the  calculated 
electron  number  density  as  a  function  of  radial  distance.  This  plot 
shows  that  a  wide  scatter  exists  in  the  electron  density  for  a  given 
radial  distance.  This  scatter  is  not  unexpected  since  the  examination 
of  individual  passes  shows  large  variations  from  orbit  to  orbit,  as 
well  as  substantial  variations  on  individual  passes.  The  scatter  plot 
also  combines  data  over  a  wide  range  of  magnetic  latitudes,  so  that 
some  of  the  scatter  could  be  due  to  variations  of  the  latitudinal 
position  of  the  spacecraft  within  the  polar  cap.  Despite  the  wide 
spread  in  the  electron  density  values  at  a  given  radial  distance  the 
electron  density  is  clearly  seen  to  decrease  with  increasing  radial 
distance,  as  would  be  expected. 

The  sharply  defined  upper  density  cutoff  evident  in  Figure  9  from 
about  3.5  to  4.6  Rg  corresponds  to  the  electron  gyrofrequency.  This 
upper  cutoff  occurs  because  the  electron  number  densities  can  be 


determined  from  the  whistler  mode  cutoff  frequency  only  when  the  elec¬ 
tron  plasma  frequency  is  lower  than  the  electron  gyrof requency.  When 
the  electron  plasma  frequency  exceeds  the  gyrof requency,  the  electron 
density  calculated  from  the  cutoff  frequency  (now  the  electron  gyrofre- 
quency)  merely  represents  a  lower  limit  on  the  actual  electron  density. 
A  statistical  tabulation  of  the  data  points  for  which  the  cutoff  fre¬ 
quency  was  greater  than  90%  of  the  electron  gyrof requency  indicated 
that  the  electron  gyrofrequency  contaminated  the  calculation  of  elec¬ 
tron  densities  in  less  than  1%  of  the  cases  below  3.5  Rg,  increasing  to 
about  3.6%  of  the  cases  for  radial  distances  above  4.25  Rg.  A  total  of 
approximately  2%  of  all  density  values  are  affected. 

To  further  explore  the  relationship  between  electron  density  and 
radial  distance  and  to  determine  the  spread  in  number  density  values 
for  a  given  altitude,  the  radial  range  was  divided  into  radial  incre¬ 
ments  of  0.25  Rg  and  the  median  value  of  the  electron  number  density 


was  computed  for  each  radial  increment.  All  cutoff  frequencies  were 
used  in  calculating  the  median  values,  even  when  the  cutoff  is  deter¬ 
mined  by  the  gyrofrequency.  Median  values  were  chosen,  rather  than 
average  values,  to  eliminate  the  effect  of  plasma  frequencies  greater 
than  the  gyrofrequency.  Since  the  medians  are  in  all  cases  well  below 
the  gyrofrequency,  median  values  are  not  affected  by  the  gyrofrequency 
cutoff.  Figure  10  is  a  log-log  plot  of  the  median  electron  number  den¬ 
sity  as  a  function  of  radial  distance.  The  error  bars  Indicate 
quartlles,  which  show  the  spread  in  25%  of  the  data  on  either  side  of 
the  median  values.  The  median  number  densities  range  from  35.2  ±  8.5 
cm t"3  at  2.1  Rg  to  0.99  ±  0.51  cm-3  at  4.66  Rg.  Variations  up  to  a 


22 


factor  of  4  occur  from  these  median  values.  As  can  be  seen  from  Figure 


V.  COMPARISON  OF  THE  DENSITY  PROFILES  WITH  THEORETICAL  MODELS 


In  Che  high  density  low-latitude  region  of  the  Earth's  magneto¬ 
sphere,  plasma  is  confined  on  closed  magnetic  field  lines  and  exists 
in  a  state  of  near-diffusive  equilibrium  [Dessler  and  Michel,  1966; 
Dungey,  1961].  In  this  region  diffusive  equilibrium  solutions  can  be 
obtained  assuming  a  Maxwellian  velocity  distribution  and  a  collision- 
dominated  plasma  [Bauer,  1969].  However,  at  high  geomagnetic  latitudes 
the  magnetic  field  lines  are  open  to  interplanetary  space.  In  this 
region  the  boundary  plasma  pressure  is  very  low  and  an  outward  flow  of 
ionospheric  ions  results  [Banks  and  Holzer,  1968,  1969b].  The  outward 
flux  of  predominantly  light  ions  has  been  called  the  polar  wind 
[Axford,  1968].  Because  the  plasma  escapes  freely  from  the  magneto¬ 
sphere  the  electron  density  drops  to  very  low  values  at  all  polar  lati¬ 
tudes.  Evidence  of  a  substantial  depletion  of  the  ionospheric  plasma 
over  the  polar  region  at  low  altitudes  has  been  given  by  Hagg  [1967] 
and  Timleck  and  Nelms  [1969]. 

The  existence  of  an  upward,  field-aligned  flux  of  protons,  pre¬ 
dicted  by  the  polar  wind  theory,  was  confirmed  by  Explorer  XXXI.  The 
ion  mass  spectrometer  on  this  satellite  detected  an  upward  proton  flux 
in  the  low  plasma  density  region  of  the  polar  cap  above  1000  km 
[Hoffman,  1969;  Hoffman  and  Dodson,  1980].  At  the  higher  altitudes  of 
the  DE-1  orbit,  the  Retarding  Ion  Mass  Spectrometer  and  the  Energetic 
Ion  Composition  Spectrometer  have  both  detected  a  similar  flux  of 
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lighter  ions  (H+  and  He+)  in  the  polar  cap  region  [Chappell  et  al., 
1982a;  Shelley  et  al.,  1982].  A  significant  flux  of  N+  ions  [Chappell 
et  al.,  1982b]  has  also  been  detected,  and  recent  measurements  have 
discovered  a  peak  in  the  flux  of  the  0+  ion  [Shelley  et  al.,  1982] 
which  appears  to  be  associated  with  periods  of  increased  magnetic 
activity  (kp  >  5)  [J.  H.  Waite,  personal  communication]. 

The  escape  of  ions  from  the  polar  cap  region  is  constrained  by  the 
continuity  equation  which,  for  a  specific  ion  species,  can  be  written 


3n 

—  +  $*(n± 


where  n*  is  the  ion  number  density,  ft*  is  the  ion  flow  velocity,  is 
the  ion  production  rate  and  %  is  the  ion  loss  rate.  Adopting  steady 
state  conditions,  the  first  term  can  be  assumed  to  be  zero.  Although 
steady  state  conditions  do  not  exist  at  all  times  in  the  polar  magneto¬ 
sphere,  steady  state  is  a  justifiable  first-order  approximation  when 
the  density  measurements  are  time-averaged  over  several  months,  as  we 
have  done.  For  the  low  neutral  particle  densities  encountered  in  the 
polar  magnetosphere  at  radial  distances  greater  than  2  Rg,  the  ion  loss 
and  production  terms  can  be  assumed  to  be  negligible  so  that  the 
continuity  equation  becomes 
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The  terms  in  parentheses  refer  to  the  ion  flux  parallel  ( n)  and  per¬ 
pendicular  (1)  to  the  magnetic  field.  The  constancy  of  the  magnetic 
field  in  the  polar  cap  requires  a  zero  divergence  of  the  perpendicular 
flow  velocity.  Although  latitudinal  variations  in  the  number  density 
could  produce  a  non-zero  divergence  in  the  perpendicular  flux,  we 
assume,  as  a  first  order  approximation,  that  this  divergence  is 
negligibly  small.  With  these  approximations,  the  continuity  equation 
for  a  simple  radial  outflow  model  requires  that 


ni  ui  A  -  constant 


where  uj  is  the  ion  flow  velocity  parallel  to  the  magnetic  field  and  A 
is  the  cross-sectional  area  of  the  flux  tube  at  a  given  radial 
distance. 

Through  the  use  of  the  conservation  of  magnetic  flux  and  a  dipole 
approximation  for  the  magnetic  field  Intensity  in  the  polar  cap  region 
(B  a  R”3),  it  can  be  shown  that  the  cross-sectional  area  of  the  ion 
flux  tube  is  proportional  to  the  cube  of  the  radial  distance.  The 
continuity  equation  for  the  steady  state  radial  outflow  of  ions  along 
magnetic  field  lines,  neglecting  ion  loss  and  production,  then  predicts 
that 
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where  A  a  R3  has  been  substituted  for  the  cross-sectional  area  of  the 


flux  tube  in  the  continuity  equation.  The  condition  of  electrical 
neutrality  in  the  polar  cap  region  requires  that 

ne  "  £  zi  ni 

where  is  the  charge  state  of  the  individual  ion  species.  Since  one 
singly-charged  ion  species  is  expected  to  dominate  at  large  radial 
distances  (H*  in  the  model  developed  by  Banks  and  Holzer),  the  electron 
density  is  controlled  by  that  species.  Assuming  therefore  that  rig  »»  n^ 
it  follows  that 


n  a  u_1  R~3 
e 

where  u  is  the  plasma  flow  velocity  along  the  magnetic  field  line.  A 
constant  plasma  flow  velocity  would  place  an  upper  limit  on  the  expo¬ 
nent  of  a  power  law  distribution  of  -3: 

rig  a  R-3  for  u  **  constant 

This  power  law  distribution  is  the  asymptotic  limit  of  the  polar  wind 
solution  derived  by  Banks  and  Holzer  [1968]. 

In  practice  one  expects  the  steady-state  flow  velocity  to  increase 
with  increasing  radial  distance  because  of  the  pressure  gradient 


between  the  ionosphere  and  the  magnetotail.  This  increase  in  the  flow 
velocity  with  Increasing  radial  distance  can  be  seen  in  the  calcula¬ 
tions  of  Banks  and  Holzer  [ 1969a, b]  and  of  Schunk  and  Watkins  [1981, 
1982].  Thus,  the  exponent  greater  than  3  in  Figure  10  is  consistent 
with  the  expected  relationships  of  radial  distance  with  flow  velocity 
and  cross-sectional  area  in  the  polar  cap  region.  It  is  reasonable  to 

— v 

expect  that  ng  a  R  where  y  >  3,  if  a  valid  power  law  relationship 
exists  between  electron  number  density  and  radial  distance.  Our  best 
data  fit  in  Figure  10,  iig  a  R~3,85  -  °*32>  indicates  that  the  plasma 
flow  velocity  increases  as  R0,85  *  0.32  (i.e.,  approximately  linearly 
with  radial  distance). 

The  newly  discovered  dominance  of  the  oxygen  ion  for  some  DE-1 
polar  passes  was  unexpected  [Shelley  et  al.,  1982].  The  classic  polar 
wind  theory  predicted  a  large  flux  of  lighter  ions  along  open  magnetic 
field  lines  at  high  latitudes  and  a  negligible  loss  of  heavy  ions. 
Calculations  showed  that  the  loss  of  light  ions  from  the  topside  polar 
ionosphere  would  result  in  a  dominance  of  0*  ions  up  to  an  altitude  of 
4000  km  to  7500  km.  Above  this  altitude  a  substantial,  though  rapidly 
diminishing,  upward  0+  flux  would  exist  to  support  the  dominant  flux 
through  charge-exchange  reactions. 

In  their  theoretical  study  of  multicomponent  plasma  transport  at 
high  geomagnetic  latitudes.  Banks  and  Holzer  [1969b]  derived  a  set  of 
solutions  to  the  standard  hydrodynamic  transport  equations  for  various 
neutral  atmosphere  models  and  a  range  of  isotropic,  equal  electron  and 
ion  temperatures.  At  an  altitude  of  7000  km  (2.18  Rg) ,  their  predicted 


electron  densities  ranged  from  approximately  80  cm-3  to  nearly  1000 
cm-3.  These  densities  exceed  our  observed  densities  at  2.25  Rg  by  a 
factor  of  5  to  50. 

Banks  and  Holzer  [1969b]  also  derived  a  set  of  predicted  ion  den¬ 
sities  for  different  neutral  atmosphere  temperatures  and  a  range  of 
equal  electron  and  ion  temperatures.  Predicted  O'**  densities  were  found 
to  vary  from  10  cm-3  to  200  cm-3  and  IT*"  densities  from  90  cm-3  to  120 
cm~3  at  7000  km.  The  altitude  of  the  transition  from  dominant  0+  at 
low  altitudes  to  dominant  H+  at  higher  altitudes  was  found  to  vary 
between  4000  km  and  7500  kmt  depending  on  the  initial  concentration  of 
hydrogen  in  the  neutral  atmosphere  model. 

When  electron  and  ion  temperatures  were  allowed  to  vary  individ¬ 
ually  and  anisotropic  temperature  distributions  were  included  in  the 
analysis  of  high-latitude  plasma  transport,  the  predicted  ion  density 
profiles  were  altered  and  the  transition  to  a  dominance  of  the  H+  ion 
was  discovered  to  occur  at  even  lower  altitudes.  In  1982,  Schunk  and 
Watkins  [1982]  developed  a  set  of  steady-state  solutions  for  the 
high-latitude  plasma  flow  in  the  altitude  range  of  1500  km  to  12,000 
km.  Their  solutions  were  derived  from  the  13-moment  system  of  trans¬ 
port  equations  which  included  continuity,  momentum,  internal  energy, 
stress  tensor  and  heat  flow  equations  for  electrons  and  oxygen  and 
hydrogen  ions.  The  solutions  were  determined  for  a  range  of  lower 
boundary  conditions  which  included  the  H+  drift  velocity,  electron  and 
H*  temperature  gradients,  ion  number  densities,  electron  and  IT*"  temper¬ 
atures  and  heat  fluxes  for  the  electrons  and  ions  at  1500  km.  Isis-2 
provided  the  measurements  of  the  hydrogen  ion  number  density,  drift 
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velocity  and  temperature  at  the  base  reference  level  and  S3-3  provided 
electron  temperature  measurements  and  a  range  of  possible  electron 
temperature  gradients  at  the  reference  level. 

Schunk  and  Watkins  [1982]  published  ion  density  predictions  both 
for  supersonic  (H  >  1)  and  subsonic  (M  <  1)  H*  outflows  where  the  Mach 
number  M  refers  to  the  ratio  of  ion  flow  velocity  to  the  thermal  velo¬ 
city.  The  13-moment  transport  equations  predict  a  significant  electron 
temperature  anisotropy  (Te^  >  Te  |)  for  the  supersonic  IT*"  outflow  and  an 
Isotropic  electron  temperature  distribution  at  this  altitude  range  for 
a  subsonic  flow.  Both  sets  of  ion  density  predictions  are  presented  in 
Figure  11  for  two  lower  boundary  electron  temperature  gradients 
(VTe  -  0. 1°K  km-1,  VTe  •  18K  km-1).  As  the  boundary  electron 
temperature  gradient  increases,  the  solutions  for  the  supersonic  H+ 
outflow  8 how  a  decrease  in  the  hydrogen  ion  density  and  a  slight 
increase  in  the  oxygen  ion  density.  Charge  neutrality  demands  that 
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and  it  can  be  seen  that  the  sum  of  the  two  highest  ion  densities  (H* 
and  0*)  predicted  by  these  steady-state  solutions  to  the  plasma  trans¬ 
port  equations  for  supersonic  H-*-  outflow  is  within  ±25%  of  the  median 
electron  density  values  observed  by  DE-1  above  2.5  Rg.  Predicted  ion 
densities  for  subsonic  H+  outflow,  on  the  other  hand,  are  much  higher 
at  2  Kg  than  our  data  would  suggest  and  the  ion  density  profile  drops 
off  with  increasing  radial  distance  much  more  slowly  than  our  mea¬ 
surements  indicate. 
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Recent  measurements  by  the  DE-1  High  Altitude  Plasma  Instrument 
have  established  the  existence  of  a  low  energy,  field-aligned  super¬ 
sonic  ion  flux  in  the  polar  cap  and  a  conic  component  of  the  ion  flux 
due  to  the  perpendicular  heating  of  the  field-aligned  component 
[Gurgiolo  and  Burch,  1982].  The  low  energy,  field-aligned  ions  consti¬ 
tute  the  classical  polar  wind  as  described  by  Axford  [1968].  These 
measurements  suggest  that  the  subsonic  H+  outflow  solutions  illustrated 
in  Figure  11  are  not  appropriate  in  the  polar  magnetosphere  at  the 
altitude  of  the  DE-1  orbit. 

Below  2.5  Rg  the  slightly  lower  theoretical  ion  densities  based  on 
the  polar  wind  model,  illustrated  in  Figure  11  for  supersonic  H+  out¬ 
flows,  are  not  unexpected  since  they  are  derived  from  the  steady-state 
solutions  to  the  plasma  transport  equations.  The  existence  of  steady- 
state  conditions  along  polar  magnetic  field  lines  is  related  to  the 
convection  of  plasma  across  the  polar  magnetosphere  and  to  the  speed  of 
the  magnetic  convective  cycle  [Banks  and  Holzer,  1969b].  Regions  of 
recently  disconnected  field  lines,  where  the  time  delay  has  not  been 
sufficient  to  achieve  a  steady  state,  can  be  expected  to  have  higher 
number  densities  than  regions  deep  in  the  polar  cap.  Thus,  the  number 
density  can  be  expected  to  vary  with  magnetic  latitudes  across  the 
polar  cap.  A  rapid  convective  speed,  which  shortens  the  time  Interval 
between  the  merging  and  reconnection  of  magnetic  field  lines,  may  be 
sufficiently  great  to  prevent  the  plasma  from  achieving  a  steady-state 
flow  in  the  polar  magnetosphere.  Under  these  circumstances,  observed 
number  densities  across  the  polar  cap  will  be  generally  higher  than  the 
number  densities  predicted  by  the  steady-state  solutions.  Variations 


in  Che  number  density  due  to  changes  in  the  magnetic  latitude  or 
disturbances  in  the  magnetic  field  on  small  time  scales  have  been 
minimised  in  this  study  by  averaging  densities  over  all  polar  magnetic 
latitudes  for  L-shell  values  greater  than  10  at  any  given  radial 
distance.  The  resulting  density  profile  is  expected  to  be  higher  than 
the  predicted  steady-state  densities. 


VI.  COMPARISON  OF  THE  HIGH-ALTITUDE  DENSITY  PROFILE  WITH 


TOPSIDE  SOUNDER  MEASUREMENTS 


Figure  12  is  a  log-log  plot  of  the  electron  number  density  as  a 
function  of  radial  distance  from  1.05  Re  to  4.7  Rg.  The  electron  den¬ 
sities  for  geocentric  distances  of  1.05  Rg  to  1.47  R£  are  median  den¬ 
sity  values  calculated  from  a  random  sampling  of  more  than  350  hlgh- 
latltude  (A  >  80°)  density  profiles  derived  from  the  Alouette  II  and 
lsis-1  topside  sounder  data.  Error  bars  Indicate  the  spread  in  elec¬ 
tron  density  values  for  this  random  sampling  of  Alouette  II  profiles 
(1966-1967)  and  Isis-l  profiles  (1969).  The  large  density  variations 
are  due  to  latitudinal  variations  in  temperature  and  ionic  composition 
[Nelms  and  Lockwood,  1967]  as  well  as  variations  in  solar  activity  and 
magnetic  field  conditions. 

An  abrupt  change  in  the  nature  of  the  Alouette/lsis  density 
profile  occurs  at  approximately  1000  km  (1.16  Rg).  This  transition 
indicates  a  change  in  the  basic  processes  of  ion  transport  in  the  lower 
ionosphere.  Below  1000  km,  conditions  of  near  diffusive  equilibrium 
combined  with  an  outward  flux  of  plasma  at  subsonic  transport  speeds 
are  expected  to  exist  [Bauer,  1969].  A  diffusive  equilibrium  model, 
based  on  the  rapid  diffusion  of  plasma  above  the  F2  peak,  predicts  an 
electron  density  distribution  of  [Angerami  and  Thomas,  1964] 


for  z  "  R-l 


and  hj  »  KT^/mig 

where  tig  is  the  electron  number  density,  R  is  the  geocentric  radial 
distance  in  Earth  radii,  k  Is  Boltzmann's  constant,  h^  is  the  scale 
height  of  the  dominant  ion  species ,  is  the  ion  temperature  and  m^  is 
the  ion  mass.  An  analysis  of  the  Alouette/lsis  density  data  below  1000 
km  yields  a  plasma  scale  height  of  337  km.  This  value  for  the  plasma 
scale  height  compares  favorably  with  the  0+  scale  height  of  282  km  for 
a  representative  ion  temperature  of  4000° K  at  an  altitude  of  1000  km. 

0+  is  the  dominant  ionospheric  ion  at  all  altitudes  below  1000  km  and 
is  the  leading  term  in  the  determination  of  the  plasma  scale  height. 

Although  a  vertical  diffusive  equilibrium  distribution  is  a  good 
approximation  to  the  Alouette/lsis  density  profiles  below  1000  km,  the 
applicability  of  the  model  is  limited  in  the  polar  cap  region  by  the 
existence  of  geomagnetic  field  lines  that  open  into  regions  of  low 
plasma  pressure  in  the  magnetotail  [Bauer,  1969].  Such  boundary  condi¬ 
tions  result  in  an  outward  plasma  flow  along  the  field  lines  at 
initially  subsonic  transport  speeds  below  1000  km  [Banks  and  Holzer, 
1969b].  In  their  solutions  of  the  hydrodynamic  transport  equations  for 
If*  and  0*,  Raitt  et  al.  [1975]  determined  that  energy  losses  due  to 
convective  plasma  motion  are  significant  in  determining  density  distri¬ 
butions  between  700  km  and  1000  km. 
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Above  1000  km,  a  weakening  of  the  charge-exchange  reactions 
between  molecules  and  Ions  and  a  decrease  In  ion-ion  friction  due  to 
rapidly  decreasing  densities  results  in  an  upward  plasma  flow  along 
magnetic  field  lines  [Banks  and  Holzer,  1969b].  This  outward  flow, 
even  for  low-energy,  field-aligned  ions,  has  been  found  to  be  super¬ 
sonic  [Gurglolo  and  Burch,  1982].  The  transition  to  supersonic  trans¬ 
port  speeds  is  theoretically  expected  to  occur  between  1000  km  and  1500 
km  [Raitt  et  al.,  1975].  However,  the  assumption  of  a  collisionless 
plasma  valid  at  the  altitude  of  the  DE-1  orbit,  is  not  valid  below  3000 
km  due  to  the  substantially  higher  densities  in  this  region.  Ion 
transport  processes,  consequently,  are  found  to  depend  in  a  complicated 
way  on  gradients  in  the  plasma  parameters  and  Coulomb  collision  fre¬ 
quencies  [Schunk  and  Watkins,  1981;  1982;  Raitt  et  al.,  1975]. 

In  the  data  gap  which  occurs  between  the  Alouette/lsis  orbits  and 
the  DE-1  orbit  (see  Figure  12),  another  transition  must  occur  in  order 
to  explain  the  change  in  the  slope  of  the  density  profiles  at  low  and 
high  altitudes.  The  observed  change  in  the  density  profile  between 
1.55  Rg  and  2.0  Rg  is  attributed  to  the  transition  from  a  collision- 
dominated  regime  below  to  a  collisionless  regime  above.  In  addition  to 
changes  in  the  basic  plasma  transport  processes,  transitions  in  the 
density  profile  may  also  be  attributed  to  changes  in  the  ionic  composi¬ 
tion  and  ion  acceleration  processes  in  the  polar  magnetosphere. 

The  topside  sounder  data  used  to  derive  the  median  density  values 
observed  in  Figure  12  are  biased  toward  regions  of  higher  densities  (> 
300  cm”3)  by  the  standard  methods  of  deriving  electron  densities  from 
topside  ionograms  [Nelms  and  Lockwood,  1967].  Electron  density  values 


35 


for  these  profiles  were  calculated  indirectly  from  the  extraordinary 
wave  reflection  frequency  and  the  gyrofrequency  and  upper  hybrid 
resonance  frequency  [Hagg,  1967;  Alouette  II  Ionospheric  Data  from 
Communications  Research  Centre  in  Ottawa]  for  ionograms  on  which  all 
traces  were  distinctly  visible.  This  indirect  method  of  deriving 
density  values  was  employed  because  the  ordinary  wave  reflection  trace 
at  the  plasma  resonance  (fp)  was  not  always  observed,  especially  in  low 
density  regions  at  high  latitudes  and  altitudes  where  the  plasma 
frequency  occasionally  fell  below  the  low  frequency  limit  of  the 
satellite's  topside  sounder. 

In  regions  of  very  low  electron  densities  and  subsequently  low 
plasma  frequencies  where  fp  «  fg,  the  upper  hybrid  resonance  frequency 
approaches  the  electron  gyrofrequency  and  the  indirect  determination  of 
electron  density  described  above  could  not  be  used.  Interference 
between  the  electron  gyrofrequency  and  the  upper  hybrid  resonance 
results  in  a  difference  or  modulation  frequency  that  can  be  determined 
from  the  ionograms  and  used  to  calculate  electron  densities  below  100 
cm-3  [Hagg,  1967].  The  use  of  a  modulation  frequency  technique  to 
determine  low  density  values  resulted  in  the  discovery  of  very  low 
electron  densities  (<  100  cm-3)  for  all  polar  latitudes  (\n  >  60°)  at 
radial  distances  of  1.24  Rg  to  1.47  Rg  for  selected  ionograms  on  which 
the  extraordinary  wave  reflection  trace  was  not  distinctly  visible 
[Nelms  and  Lockwood,  1967;  Hagg,  1967].  The  inclusion  of  these  very 
low  densities  on  Figure  12  would  lower  the  Alouette/  Isis  density  curve 
by  an  order  of  magnitude  above  1.2  Rg.  This  would  reduce  the  discre¬ 
pancy  between  the  low-altitude  density  studies  and  the  observed  DE-1 


density  values,  since  the  whistler  mode  frequency  cutoff  method  of 
determining  electron  densities  is  valid  for  low  densities  (see  Figure 
10).  A  more  complete  analysis  of  available  topside  ionograms  would  be 
necessary  in  order  to  determine  an  accurate  slope  of  the  electron 
density  profile  between  1.24  Rg  and  1.47  Rg. 
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VII.  DISCUSSION  AND  SUMMARY 


Electron  densities  at  2.0  Rg  to  4.66  Rg  have  been  derived  from  the 
plasma  frequency  cutoff  of  the  whistler  mode  radiation  in  the  polar  cap 
using  electric  field  spectrum  measurements  from  the  Plasma  Wave  Instru¬ 
ment  on  the  DE-1  spacecraft.  A  log-log  plot  of  the  median  density 
values  as  a  function  of  radial  distance  indicates  a  power  law  distribu¬ 
tion  of  electron  densities,  varying  as  R-3*85  ±  0.32.  This  power  law 
distribution  is  consistent  with  a  steady  state,  radial  outflow  of  iono¬ 
spheric  ions  and  electrons  along  polar  magnetic  field  lines.  It 
implies  a  nearly  linear  increase  in  plasma  flow  velocity  with  increas¬ 
ing  radial  distance  at  radial  distances  above  2  Rg. 

A  comparison  of  observed  electron  densities  with  theoretical  ion 
and  electron  densities  based  on  the  classic  polar  wind  model  for  super¬ 
sonic  ion  flow  velocities  yields  consistent  results  at  2.5  Rg  to  2.8 
Rg.  The  theoretical  polar  wind  densities  below  2.5  Rg  are  slightly 
lower  than  the  observed  median  densities .  This  result  is  not  unex¬ 
pected  since  conditions  of  steady-state  plasma  flow  are  not  expected  to 
exist  across  the  polar  cap  at  all  latitudes  and  at  all  times. 

A  comparison  of  observed  electron  densities  with  low-altitude 
density  profiles  from  the  Alouette/Isis  topside  sounder  data  Indicates 
several  evolutions  in  the  radial  dependence  of  electron  distributions 
in  the  polar  cap.  A  state  of  near  diffusive  equilibrium  combined  with 
an  outward  flux  of  plasma  at  subsonic  transport  speeds  below  1000  km 
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changes  to  a  collision-dominated  outward  flux  of  plasma  with  a  transi¬ 
tion  to  supersonic  velocities  above  1000  km.  This  transition  is 
clearly  evident  in  the  low-altitude  density  profile  at  approximately 
1.16  Rg  (see  Figure  12).  A  second  transition  is  also  indicated  in  the 
density  data  gap  between  1.55  Rg  and  2.0  Rg  to  account  for  the  changing 
nature  of  the  density  profiles  at  high  and  low  altitudes  in  the  polar 
cap.  This  second  transition  in  the  density  profile  is  believed  to  be 
due  to  a  change  from  a  collision-dominated  radial  outflow  distribution 
model  below  1.5  Rg  to  a  collisionless  power  law  distribution  above  2.0 
Rg.  Other  factors  affecting  the  change  in  the  radial  dependence  of  the 
electron  number  density  Include  changes  in  the  ionic  composition  of  the 
plasma  as  well  as  changes  in  the  plasma  flow  velocity  and  plasma 
temperature. 

Large  variations  seen  in  low-altitude  electron  densities  in 
earlier  studies  persist  at  the  higher  altitudes  of  the  DE-1  orbit.  The 
spread  in  the  density  data  is  the  result  of  averaging  over  all  polar 
magnetic  latitudes  and  nearly  80  polar  passes  for  the  five-month  time 
interval  of  this  study.  Such  a  survey  would  conceal  all  density  irreg¬ 
ularities  on  a  short  time  scale  and  would  even  suppress  irregularities 
of  large  spatial  dimensions  which  are  not  consistent  across  the  polar 
cap  at  a  constant  radial  distance.  Such  irregularities  in  the  density 
profile  at  low  altitudes  were  discovered  by  the  early  topside  sounder 
satellites.  Alouette  1,  in  a  polar  orbit  at  a  constant  radial  distance 
of  1.16  Rg,  provided  evidence  of  latitudinal,  diurnal  and  seasonal 
variations  in  the  electron  density  profile  at  high  geomagnetic  lati¬ 
tudes  [Nelms,  1966;  Thomas  et  al.,  1966].  Superimposed  on  the  mono- 
tonic  decrease  in  number  density  with  increasing  latitude  were 


prominent  peaks  and  troughs  In  the  density  distribution  on  spatial 
scales  of  several  hundred  km  [Thomas  et  al.,  1966].  In  addition,  many 
small  scale  irregularities  in  space  and  time  contributed  to  the  spread 
in  the  density  data.  Some  large  scale  irregularities  persisted  over 
continuous  polar  passes  in  narrow  bands  of  invariant  latitude  [Sato  and 
Colin,  1969;  Chan  and  Colin,  1969]..  In  general  polar  cap  density  dis¬ 
tributions  were  found  to  be  very  irregular  and  to  sensitively  depend  on 
magnetic  activity,  seasonal  variations  and  local  time. 

Similar  irregularities  in  the  density  profile  are  expected  to  per¬ 
sist  at  DE-1  altitudes.  The  previously  cited  examples  of  density  pro¬ 
file  variability  between  successive  polar  passes  (see  Figures  4-8) 
are  evidence  of  large  density  fluctuations  in  time  periods  of  less  than 
8  hours.  Such  a  large  transformation  in  the  density  profile  for  the 
same  region  of  the  polar  magnetosphere  is  probably  due  to  abrupt 
changes  in  the  magnetic  field  activity.  Sufficient  data  to  examine 
seasonal  variations  in  the  density  profile  at  DE-1  altitudes  are 
not  yet  available.  Statistically  significant  density  data  for  seasonal 
comparisons  at  the  same  radial  distances  in  the  polar  magnetosphere 
cannot  be  obtained  until  the  DE-1  orbit  has  precessed  around  to  the 
northern  polar  cap  region  in  the  spring  of  1985.  But  sufficient  den¬ 
sity  data  already  exist  to  investigate  latitudinal  density  variations 
in  the  polar  cap  region  and  to  compare  these  results  with  low-altitude 
studies  of  latitudinal  density  variations  in  order  to  determine  the 
radial  extent  of  persistent  maxima  and  minima  in  the  electron  density 
distributions.  A  further  comparison  of  magnetic  activity  and  the  den¬ 
sity  profile  as  a  function  of  time  is  planned  in  order  to  determine  the 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


The  change  in  the  orientation  of  the  DE-1  orbit  from 
September,  1981  to  February,  1982.  All  density  data 
presented  in  this  study  for  the  five-month  interval, 
indicated  by  the  cross  shading,  are  from  high- 
latitude  passes  over  the  polar  cap  region. 

A  representative  spectrogram  of  electric  field 
amplitude  measurements  illustrating  the  various 
plasma  wave  emissions  found  in  the  polar  cap  region 
near  the  electron  gyrofrequency  and  plasma 
frequency. 

An  electric  field  spectrogram  of  a  polar  cusp 
crossing  in  early  October.  Irregularities  in  the 
density  profile  are  indicated  by  sudden  Increases  in 
fp  at  0752  and  at  0845-0850  as  the  satellite  nears 
the  pole. 

A  nightside  crossing  of  the  auroral  field  lines  in 
mid-February  illustrating  a  depletion  in  electron 
densities  at  “  50°  [Calvert,  1981]  and  a  smooth 
density  profile  inside  the  polar  cap  region. 


Figure  5  A  continuation  of  the  same  polar  pass  illustrated  in 

Figure  4.  The  density  profile  is  smoothly  varying 
over  the  entire  polar  cap  region. 

Figure  6  A  spectrogram  illustrating  the  time  variability  of 

the  density  profile.  This  pass  preceded  the  pass 
illustrated  in  Figures  4  and  S  by  just  6  hours  for 
the  same  region  of  the  polar  magnetosphere.  The  den¬ 
sity  profile  is  highly  variable  on  very  short  spatial 
scales. 

Figure  7  Another  nightside  crossing  of  the  auroral  field  lines 

in  early  February.  Changes  in  the  cutoff  frequency 
for  this  pass  indicate  small  variations  in  the  number 
density  by  a  factor  of  2  or  less. 

Figure  8  This  pass  follows  the  pass  illustrated  in  Figure  7  by 

7.5  hours  for  the  region  of  the  polar  magnetosphere 
just  poleward  of  the  nightside  auroral  zone.  The 
density  profile  is  highly  irregular  with  fluctuations 
up  to  an  order  of  magnitude  on  spatial  scales  of  less 
than  700  km. 

Figure  9  A  scatter  plot  of  electron  density  as  a  function  of 

radial  distance.  The  large  scatter  in  the  data  is 
mainly  due  to  temporal  variations  in  the  electron 
density.  The  sharp  upper  cutoff  evident  from  about 
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Figure  10 


Figure  11 


Figure  12 


3.5  to  4.6  Rg  Is  caused  by  the  whistler  mode  propaga¬ 
tion  cutoff  at  the  electron  gyrofrequency. 

A  log-log  plot  of  median  number  densities  as  a  func¬ 
tion  of  radial  distance.  The  best  line  fit  through 
these  points  indicates  a  power  law  distribution  for 
electron  densities  in  the  polar  cap:  o  R“3,85. 

Predicted  ion  densities  below  2.8  Rg  based  on  solu¬ 
tions  to  the  13-moment  transport  equations  [Schunk 
and  Watkins,  1982].  The  subsonic  outflow  solutions 
correspond  to  higher  H+  and  0^  densities.  These  pre¬ 
dicted  densities  exceed  the  DE-1  measured  densities 
at  2  Rg  by  more  than  an  order  of  magnitude.  Super¬ 
sonic  outflow  solutions  correspond  to  lower  IT1*  and  0+ 
densities.  An  increase  in  the  boundary  electron 
temperature  gradient  (from  0. 1°K  km-1  to  1°K  km-1) 
corresponds  to  a  slight  increase  in  the  predicted  CT*" 
densities  and  a  larger  decrease  in  the  predicted 
densities.  Both  sets  of  predicted  ion  densities  for 
the  supersonic  outflow  solutions  are  within  25%  of 
the  median  density  values  observed  by  DE-1. 

A  log-log  plot  comparing  low-altitude  median  electron 
densities  obtained  from  the  topside  sounder  data  of 
Alouette  II  and  Isis-1  and  high-altitude  median 
electron  densities  obtained  from  the  PWI  on  DE-1. 
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